Wetland Flowpaths Mediate Nitrogen and Phosphorus Concentrations across the Upper Mississippi River Basin.

Eutrophication, harmful algal blooms, and human health impacts are critical environmental challenges resulting from excess nitrogen and phosphorus in surface waters. Yet we have limited information regarding how wetland characteristics mediate water quality across watershed scales. We developed a large, novel set of spatial variables characterizing hydrological flowpaths from wetlands to streams, that is, "wetland hydrological transport variables," to explore how wetlands statistically explain the variability in total nitrogen (TN) and total phosphorus (TP) concentrations across the Upper Mississippi River Basin (UMRB) in the United States. We found that wetland flowpath variables improved landscape-to-aquatic nutrient multilinear regression models (from R2 = 0.89 to 0.91 for TN; R2 = 0.53 to 0.84 for TP) and provided insights into potential processes governing how wetlands influence watershed-scale TN and TP concentrations. Specifically, flowpath variables describing flow-attenuating environments, for example, subsurface transport compared to overland flowpaths, were related to lower TN and TP concentrations. Frequent hydrological connections from wetlands to streams were also linked to low TP concentrations, which likely suggests a nutrient source limitation in some areas of the UMRB. Consideration of wetland flowpaths could inform management and conservation activities designed to reduce nutrient export to downstream waters.

Random Forest models to estimate bankfull and low flow channel widths and depths across the conterminous United States.

Channel dimensions (width and depth) at varying flows influence a host of instream ecological processes, as well as habitat and biotic features; they are a major consideration in stream habitat restoration and instream flow assessments. Models of widths and depths are often used to assess climate change vulnerability, develop endangered species recovery plans, and model water quality. However, development and application of such models require specific skillsets and resources. To facilitate acquisition of such estimates, we created a dataset of modeled channel dimensions for perennial stream segments across the conterminous U.S. We used random forest models to predict wetted width, thalweg depth, bankfull width, and bankfull depth from several thousand field measurements of the National Rivers and Streams Assessment. Observed channel widths varied from <5 m to >2000 m and depths varied from <2 m to >125 m. Metrics of watershed area, runoff, slope, land use, and more were used as model predictors. The models had high pseudo R-squared values (0.70 to 0.91) and median absolute errors within ±6% to ±21% of the interquartile range of measured values across ten stream orders. Predicted channel dimensions can be joined to 1.1 million stream segments of the 1:100K resolution National Hydrography Dataset Plus (version 2.1). These predictions, combined with a rapidly growing body of nationally available data, will further enhance our ability to study and protect aquatic resources.

National hydrologic connectivity classification links wetlands with stream water quality.

Wetland hydrologic connections to downstream waters influence stream water quality. However, no systematic approach for characterizing this connectivity exists. Here using physical principles, we categorized conterminous US freshwater wetlands into four hydrologic connectivity classes based on stream contact and flowpath depth to the nearest stream: riparian, non-riparian shallow, non-riparian mid-depth and non-riparian deep. These classes were heterogeneously distributed over the conterminous United States; for example, riparian dominated the south-eastern and Gulf coasts, while non-riparian deep dominated the Upper Midwest and High Plains. Analysis of a national stream dataset indicated acidification and organic matter brownification increased with connectivity. Eutrophication and sedimentation decreased with wetland area but did not respond to connectivity. This classification advances our mechanistic understanding of wetland influences on water quality nationally and could be applied globally.

Estimating biotic integrity to capture existence value of freshwater ecosystems.

The US Environmental Protection Agency (EPA) uses a water quality index (WQI) to estimate benefits of proposed Clean Water Act regulations. The WQI is relevant to human use value, such as recreation, but may not fully capture aspects of nonuse value, such as existence value. Here, we identify an index of biological integrity to supplement the WQI in a forthcoming national stated preference survey that seeks to capture existence value of streams and lakes more accurately within the conterminous United States (CONUS). We used literature and focus group research to evaluate aquatic indices regularly reported by the EPA's National Aquatic Resource Surveys. We chose an index that quantifies loss in biodiversity as the observed-to-expected (O/E) ratio of taxonomic composition because focus group participants easily understood its meaning and the environmental changes that would result in incremental improvements. However, available datasets of this index do not provide the spatial coverage to account for how conditions near survey respondents affect their willingness to pay for its improvement. Therefore, we modeled and interpolated the values of this index from sampled sites to 1.1 million stream segments and 297,071 lakes across the CONUS to provide the required coverage. The models explained 13 to 36% of the variation in O/E scores and demonstrate how modeling can provide data at the required density for benefits estimation. We close by discussing future work to improve performance of the models and to link biological condition with water quality and habitat models that will allow us to forecast changes resulting from regulatory options.

A multiscale landscape approach for prioritizing river and stream protection and restoration actions.

River and stream conservation programs have historically focused on a single spatial scale, for example, a watershed or stream site. Recently, the use of landscape information (e.g., land use and land cover) at multiple spatial scales and over large spatial extents has highlighted the importance of incorporating a landscape perspective into stream protection and restoration activities. Previously, we developed a novel framework that links information about watershed-, catchment-, and reach-scale integrity with stream biological condition using scatterplots and a landscape integrity map. Here we examined an application of this approach for streams in urban and other settings in King County, Washington State, United States, where we related stream macroinvertebrate condition to two indices of landscape integrity, the US Environmental Protection Agency's (USEPA) nationally available Index of Watershed Integrity (IWI) and Index of Catchment Integrity (ICI). We generated a scatterplot of IWI versus ICI for sample sites, where points represented site macroinvertebrate condition from poor to good. The same data were also visualized as a landscape integrity map that displayed catchments of King County according to the level of watershed and catchment integrity (high or low IWI/ICI). Almost three-quarters of poor-condition sites were associated with high-integrity watersheds and catchments (i.e., underperforming sites), which suggested that either one or both national indicators were insufficient for this area, and that sites underperformed because of local-scale factors. In response, we used a catchment-scale indicator related to forest condition (PctForestCat) after examining several GIS-based dispersal indicators from the National Hydrography Dataset and other candidates from the USEPA's StreamCat dataset. We then compared the results of the scatterplots and maps based on the current and original analyses and found that many of the sites previously classified as underperforming now performed as expected, that is, they were poor-condition sites in poor-condition catchments. This analysis demonstrates how results based on a national dataset can be improved by developing an alternative that represents regionally important stressors. The methods used to develop an effective landscape indicator based on StreamCat datasets, and the utility of the multiscale approach, could provide important tools for prioritizing, optimizing, and communicating stream conservation actions.

Identifying lakes at risk of toxic cyanobacterial blooms using satellite imagery and field surveys across the United States.

Harmful algal blooms caused by cyanobacteria are a threat to global water resources and human health. Satellite remote sensing has vastly expanded spatial and temporal data on lake cyanobacteria, yet there is still acute need for tools that identify which waterbodies are at-risk for toxic cyanobacterial blooms. Algal toxins cannot be directly detected through imagery but monitoring toxins associated with cyanobacterial blooms is critical for assessing risk to the environment, animals, and people. The objective of this study is to address this need by developing an approach relating satellite imagery on cyanobacteria with field surveys to model the risk of toxic blooms among lakes. The Medium Resolution Imaging Spectrometer (MERIS) and United States (US) National Lakes Assessments are leveraged to model the probability among lakes of exceeding lower and higher demonstration thresholds for microcystin toxin, cyanobacteria, and chlorophyll a. By leveraging the large spatial variation among lakes using two national-scale data sources, rather than focusing on temporal variability, this approach avoids many of the previous challenges in relating satellite imagery to cyanotoxins. For every satellite-derived lake-level Cyanobacteria Index (CI_cyano) increase of 0.01 CI_cyano/km2, the odds of exceeding six bloom thresholds increased by 23-54 %. When the models were applied to the 2192 satellite monitored lakes in the US, the number of lakes identified with ≥75 % probability of exceeding the thresholds included as many as 335 lakes for the lower thresholds and 70 lakes for the higher thresholds, respectively. For microcystin, the models identified 162 and 70 lakes with ≥75 % probability of exceeding the lower (0.2 µg/L) and higher (1.0 µg/L) thresholds, respectively. This approach represents a critical advancement in using satellite imagery and field data to identify lakes at risk for developing toxic cyanobacteria blooms. Such models can help translate satellite data to aid water quality monitoring and management.

To burn or not to burn: An empirical assessment of the impacts of wildfires and prescribed fires on trace element concentrations in Western US streams.

The use of low-severity prescribed fires has been increasingly promoted to reduce the impacts from high-severity wildfires and maintain ecosystem resilience. However, the effects of prescribed fires on water quality have rarely been evaluated relative to the effects of wildfires. In this study, we assessed the effects of 54 wildfires and 11 prescribed fires on trace element (arsenic, selenium, and cadmium) concentrations of streams draining burned watersheds in the western US. To obtain results independent of the choice of method, we employed three independent analytical approaches to evaluate fire effects on water quality for the first three post-fire years. In general, we observed significant increases in trace element concentrations in streams burned by large, high-severity wildfires, despite substantial variability across sites. Comparatively, we did not observe increases in the spring mean concentration of arsenic, selenium, and cadmium in watersheds burned by prescribed fires. Our analysis indicated that the post-fire trace element response in streams was primarily influenced by burn area, burn severity, post-fire weather, surface lithology, watershed physiography, and land cover. This study's results demonstrate that prescribed burns could lessen the post-fire trace element loads in downstream waters if prescribed fires reduce subsequent high severity fires in the landscape.

Variable wildfire impacts on the seasonal water temperatures of western US streams: A retrospective study.

Recent increases in the burn area and severity of wildfires in the western US have raised concerns about the impact on stream water temperature-a key determinant of cold-water fish habitats. However, the effect on seasonal water temperatures of concern, including winter and summer, are not fully understood. In this study, we assessed the impact of wildfire burns at Boulder Creek (Oregon), Elk Creek (Oregon), and Gibbon River (Wyoming) watersheds on the downstream winter and summer water temperatures for the first three post-fire years. To obtain results independent of the choice of the analytical method, we evaluated the consequence of each burn using three different statistical approaches that utilize local water temperature data. Our results from the three approaches indicated that the response of water temperatures to wildfire burns varied across seasons and sites. Wildfire burns were associated with a median increase of up to 0.56°C (Standard Error; S.E. < 0.23°C) in the summer mean water temperatures (MWT) and 62 degree-day Celsius (DDC; S.E. < 20.7 DDC) in the summer accumulated degree days (ADD) for the three subsequent years across studied stream sites. Interestingly, these burns also corresponded to a median decrease of up to 0.49°C (S.E. < 0.45°C) in the winter MWT and 39 DDC (S.E. < 40.5 DDC) in the winter ADD for the same period across sites. Wildfire effects on the downstream water temperatures diminished with increasing site distance from the burn perimeter. Our analyses demonstrated that analytical methods that utilize local watershed data could be applied to evaluate fire effects on downstream water temperatures.

Geospatial Patterns of Antimicrobial Resistance Genes in the US EPA National Rivers and Streams Assessment Survey.

Antimicrobial resistance (AR) is a serious global problem due to the overuse of antimicrobials in human, animal, and agriculture sectors. There is intense research to control the dissemination of AR, but little is known regarding the environmental drivers influencing its spread. Although AR genes (ARGs) are detected in many different environments, the risk associated with the spread of these genes to microbial pathogens is unknown. Recreational microbial exposure risks are likely to be greater in water bodies receiving discharge from human and animal waste in comparison to less disturbed aquatic environments. Given this scenario, research practitioners are encouraged to consider an ecological context to assess the effect of environmental ARGs on public health. Here, we use a stratified, probabilistic survey of nearly 2000 sites to determine national patterns of the anthropogenic indicator class I integron Integrase gene (intI1) and several ARGs in 1.2 million kilometers of United States (US) rivers and streams. Gene concentrations were greater in eastern than in western regions and in rivers and streams in poor condition. These first of their kind findings on the national distribution of intI1 and ARGs provide new information to aid risk assessment and implement mitigation strategies to protect public health.

Vulnerable Waters are Essential to Watershed Resilience.

Watershed resilience is the ability of a watershed to maintain its characteristic system state while concurrently resisting, adapting to, and reorganizing after hydrological (for example, drought, flooding) or biogeochemical (for example, excessive nutrient) disturbances. Vulnerable waters include non-floodplain wetlands and headwater streams, abundant watershed components representing the most distal extent of the freshwater aquatic network. Vulnerable waters are hydrologically dynamic and biogeochemically reactive aquatic systems, storing, processing, and releasing water and entrained (that is, dissolved and particulate) materials along expanding and contracting aquatic networks. The hydrological and biogeochemical functions emerging from these processes affect the magnitude, frequency, timing, duration, storage, and rate of change of material and energy fluxes among watershed components and to downstream waters, thereby maintaining watershed states and imparting watershed resilience. We present here a conceptual framework for understanding how vulnerable waters confer watershed resilience. We demonstrate how individual and cumulative vulnerable-water modifications (for example, reduced extent, altered connectivity) affect watershed-scale hydrological and biogeochemical disturbance response and recovery, which decreases watershed resilience and can trigger transitions across thresholds to alternative watershed states (for example, states conducive to increased flood frequency or nutrient concentrations). We subsequently describe how resilient watersheds require spatial heterogeneity and temporal variability in hydrological and biogeochemical interactions between terrestrial systems and down-gradient waters, which necessitates attention to the conservation and restoration of vulnerable waters and their downstream connectivity gradients. To conclude, we provide actionable principles for resilient watersheds and articulate research needs to further watershed resilience science and vulnerable-water management.

Wildfires can increase regulated nitrate, arsenic, and disinfection byproduct violations and concentrations in public drinking water supplies.

Wildfires are a concern for water quality in the United States, particularly in the wildland-urban interface of populous areas. Wildfires combust vegetation and surface soil organic matter, reduce plant nutrient uptake, and can alter the composition of runoff and receiving waters. At the wildland-urban interface, fires can also introduce contaminants from the combustion of man-made structures. We examine post-wildfire effects on drinking water quality by evaluating concentrations and maximum contaminant level (MCL) violations of selected contaminants regulated in the U.S. at public drinking water systems (PWSs) located downstream from wildfire events. Among contaminants regulated under the U.S. Safe Drinking Water Act, nitrate, arsenic, disinfection byproducts, and volatile organic compounds (VOCs) were analyzed in watersheds that experienced major wildfires. Surface water sourced drinking water (SWDW) nitrate violations increased by an average of 0.56 violations per PWS and concentrations increased by 0.044 mg-N/L post-wildfire. Groundwater sourced drinking water (GWDW) nitrate violations increased by 0.069 violations per PWS and concentrations increased by 0.12 mg-N/L post-wildfire. SWDW total trihalomethane (TTHM) violations increased by 0.58 violations per PWS and concentrations increased by 10.4 µg/L. SWDW total haloacetic acid (HAA5) violations increased by 0.82 violations per PWS and concentrations increased by 8.5 µg/L. Arsenic violations increased by 1.08 violations per PWS and concentrations increased by 0.92 µg/L. There was no significant effect of wildfires on average VOC violations. Nitrate violations increased in 75% of SWDW sites and 34% of GWDW sites post-wildfire, while about 71% and 50% of SWDW sites showed an increase in TTHM and HAA5 violations. Violations also increased for 35% of arsenic and 44% of VOC sites post-wildfire. These findings support the need for increased awareness about the impact of wildfires on drinking water treatment to help PWS operators adapt to the consequences of wildfires on source water quality, particularly in wildfire-prone regions.

Parsing Weather Variability and Wildfire Effects on the Post-Fire Changes in Daily Stream Flows : A Quantile-Based Statistical Approach and its Application.

Determining wildland fire impacts on streamflow can be problematic as the hydrology in burned watersheds is influenced by post-fire weather conditions. This study presents a quantile-based analytical framework for assessing fire impacts on low and peak daily flow magnitudes, while accounting for post-fire weather influences. This framework entails (a) the bootstrap method to compute the relative change in the post-fire annual flow and weather statistics, (b) Double Mass analysis to detect if post-fire baseflow and quick flow yield ratios are significantly altered, and (c) a quantile regression method to parse fire effects on flow at a specific quantile. We illustrate the applicability of this analytical framework using 44 western US streams with at least 5% of their watershed area burned. Results indicate that large, high-severity burns in upland watersheds can raise the streamflow magnitude at the 0.05 th and 0.95 th quantiles for at least the five post-fire years. Quantile regression results show that the median fire-related increase in flow for the five post-fire years can be up to 5000% (Standard Error; SE < 2%) at the 0.05 th quantile and 161% (SE < 10%) at the 0.95 th quantile. The fire-related increase in flow was often pronounced at the 0.05 th quantile for streams in the Pacific Northwest and California regions. The difference in fire effects on flow (at both quantiles) across streams was related to post-fire weather, pyrology, physiography, and land cover. The proposed analytical framework can be useful for detecting and quantifying fire effects on the low and peak stream flows in burned watersheds without overlapping disturbances.

Using hydrologic landscape classification and climatic time series to assess hydrologic vulnerability of the western U.S. to climate.

We apply the hydrologic landscape (HL) concept to assess the hydrologic vulnerability of the western United States (U.S.) to projected climate conditions. Our goal is to understand the potential impacts of hydrologic vulnerability for stakeholder-defined interests across large geographic areas. The basic assumption of the HL approach is that catchments that share similar physical and climatic characteristics are expected to have similar hydrologic characteristics. We use the hydrologic landscape vulnerability approach (HLVA) to map the HLVA index (an assessment of climate vulnerability) by integrating hydrologic landscapes into a retrospective analysis of historical data to assess variability in future climate projections and hydrology, which includes temperature, precipitation, potential evapotranspiration, snow accumulation, climatic moisture, surplus water, and seasonality of water surplus. Projections that are beyond 2 standard deviations of the historical decadal average contribute to the HLVA index for each metric. Separating vulnerability into these seven separate metrics allows stakeholders and/or water resource managers to have a more specific understanding of the potential impacts of future conditions. We also apply this approach to examine case studies. The case studies (Mt. Hood, Willamette Valley, and Napa-Sonoma Valley) are important to the ski and wine industries and illustrate how our approach might be used by specific stakeholders. The resulting vulnerability maps show that temperature and potential evapotranspiration are consistently projected to have high vulnerability indices for the western U.S. Precipitation vulnerability is not as spatially uniform as temperature. The highest-elevation areas with snow are projected to experience significant changes in snow accumulation. The seasonality vulnerability map shows that specific mountainous areas in the west are most prone to changes in seasonality, whereas many transitional terrains are moderately susceptible. This paper illustrates how HL and the HLVA can help assess climatic and hydrologic vulnerability across large spatial scales. By combining the HL concept and HLVA, resource managers could consider future climate conditions in their decisions about managing important economic and conservation resources.

Applying the index of watershed integrity to the Matanuska-Susitna basin.

The Matanuska-Susitna Borough is the fastest growing region in the State of Alaska and is impacted by a number of human activities. We conducted a multiscale assessment of the stressors facing the borough by developing and mapping the Index of Watershed Integrity (IWI) and Index of Catchment Integrity (the latter considers stressors in areas surrounding individual stream segments exclusive of upstream areas). The assessment coincided with the borough's stormwater management planning. We adapted the list of anthropogenic stressors used in the original conterminous United States IWI application to reflect the borough's geography, human activity, and data availability. This analysis also represents an early application of the NHDPlus High Resolution geospatial framework and the first use of the framework in an IWI study. We also explored how remediation of one important stressor, culverts, could impact watershed integrity at the catchment and watershed scales. Overall, we found that the integrity scores for the Matanuska-Susitna basin were high compared to the conterminous United States. Low integrity scores did occur in the rapidly developing Wasilla-Palmer core area. We also found that culvert remediation had a larger proportional impact in catchments with fewer stressors.

The use of multiscale stressors with biological condition assessments: A framework to advance the assessment and management of streams.

Incorporating information on landscape condition (or integrity) across multiple spatial scales and over large spatial extents in biological assessments may allow for a more integrated measure of stream biological condition and better management of streams. However, these systems are often assessed and managed at an individual scale (e.g., a single watershed) without a larger regional multiscale context. In this paper, our goals were: (1) To develop a conceptual framework that could combine stream biological condition to abiotic landscape integrity (or, conversely, stressor) data at three spatial scales: watershed, catchment and stream-reach scale, to enable more targeted management actions. Measures of landscape integrity and stressors are negatively related, i.e., integrity on a 0-1 scale is equal or equivalent to stressors on a 1-0 scale. (2) To develop the framework in such a way that allows operational flexibility, whereby different indicators can be used to represent biological condition, and landscape integrity (or stressors) at various scales. (3) To provide different examples of the framework's use to demonstrate the flexibility of its application and relevance to management. Examples include stream biological assessments from different regions and states across the U.S. for fish, macroinvertebrates and diatoms using a variety of assessment tools (e.g., the Biological Condition Gradient (BCG), and an Index of Biotic Integrity (IBI)). Landscape integrity indicators comprise U.S. EPA's nationally available Index of Watershed Integrity (IWI) and Index of Catchment Integrity (ICI), and state and regional derived watershed and stream-reach scale integrity indicators. Scatterplots and a landscape integrity map were used to relate samples of stream condition classes (e.g., good, fair, poor) to watershed, catchment and stream-reach scale integrity. This framework and approach could provide a powerful tool for prioritizing, targeting, and communicating management actions to protect and restore stream habitats, and for informing the spatial extent at which management is applied.

Patterns and predictions of drinking water nitrate violations across the conterminous United States.

Excess nitrate in drinking water is a human health concern, especially for young children. Public drinking water systems in violation of the 10 mg nitrate-N/L maximum contaminant level (MCL) must be reported in EPA's Safe Drinking Water Information System (SDWIS). We used SDWIS data with random forest modeling to examine the drivers of nitrate violations across the conterminous U.S. and to predict where public water systems are at risk of exceeding the nitrate MCL. As explanatory variables, we used land cover, nitrogen inputs, soil/hydrogeology, and climate variables. While we looked at the role of nitrate treatment in separate analyses, we did not include treatment as a factor in the final models, due to incomplete information in SDWIS. For groundwater (GW) systems, a classification model correctly classified 79% of catchments in violation and a regression model explained 43% of the variation in nitrate concentrations above the MCL. The most important variables in the GW classification model were % cropland, agricultural drainage, irrigation-to-precipitation ratio, nitrogen surplus, and surplus precipitation. Regions predicted to have risk for nitrate violations in GW were the Central California Valley, parts of Washington, Idaho, the Great Plains, Piedmont of Pennsylvania and Coastal Plains of Delaware, and regions of Wisconsin, Iowa, and Minnesota. For surface water (SW) systems, a classification model correctly classified 90% of catchments and a regression model explained 52% of the variation in nitrate concentration. The variables most important for the SW classification model were largely hydroclimatic variables including surplus precipitation, irrigation-to-precipitation ratio, and % shrubland. Areas at greatest risk for SW nitrate violations were generally in the non-mountainous west and southwest. Identifying the areas with possible risk for future violations and potential drivers of nitrate violations across U.S. can inform decisions on how source water protection and other management options could best protect drinking water.

Adapting the Index of Watershed Integrity for Watershed Managers in the Western Balkans Region.

Sustainable development supports watershed processes and functions. To aid the sustainable development of the western Balkans' transboundary river and lake basins, the Regional Environmental Center for Central and Eastern Europe and the US Environmental Protection Agency (EPA) adapted the EPA's Index of Watershed Integrity (IWI) following the devasting 2014 floods in Albania, Bosnia and Herzegovina, Kosovo, North Macedonia, Montenegro, and Serbia. The IWI evaluates six watershed functions based on a suite of anthropogenic stressors (e.g., impervious surfaces, reservoirs). A key feature of the IWI is its ability to accumulate the impact of upstream activities of any specific location in a river network. A novel feature of the IWI, compared with other watershed assessment tools, is its capacity to provide actionable information at the local scale. IWI scores-ranging from 0 (low integrity) to 1 (high integrity)-calculated for the 1084 catchments of the study area indicated highest integrity in the Alpine geographic region (mean = 0.55, standard deviation (SD) = 0.11) and intermediate to lowest integrity within the Mediterranean (mean = 0.49, SD = 0.12) and Continental (mean = 0.40, SD = 0.10) geographic regions. The IWI results are presented hierarchically for data analysts (stressor, functional component, Index of Catchment Integrity and IWI), ecologists (stream/catchment, watershed, basin), and managers (local, national, international). We provide real-world examples for managers, and suggestions for improving the assessment.

Watershed integrity and associations with gastrointestinal illness in the United States.

Gastrointestinal (GI) illnesses are associated with various environmental factors, such as water quality, stormwater runoff, agricultural runoff, sewer overflows, and wastewater treatment plant effluents. However, rather than assessing an individual factor alone, two indices incorporating a combination of ecological and environmental stressors were created to represent (1) overall watershed integrity, Index of Watershed Integrity (IWI) and (2) catchment integrity, Index of Catchment Integrity (ICI). These indices could provide a more comprehensive understanding of how watershed/catchment integrity potentially impact the rates of GI illness, compared to assessing an individual stressor alone. We utilized the IWI and ICI, as well as agricultural and urban land uses, to assess associations at the county level with the rates of GI illness in a population of adults over 65 years of age. Our findings demonstrated that both watershed and catchment integrity are associated with reduced hospitalizations for any GI outcomes, though association varied by urbanicity. We believe that improved versions of the IWI and ICI may potentially be useful indicators for public health analyses in other circumstances, particularly when considering rural areas or to capture the complex stressors impacting the ecological health of a watershed.

Arsenic Drinking Water Violations Decreased across the United States Following Revision of the Maximum Contaminant Level.

Arsenic poses a threat to public health due to widespread environmental prevalence and known carcinogenic effects. In 2001, the US EPA published the Final Arsenic Rule (FAR) for public drinking water, reducing the maximum contaminant level (MCL) from 50 to 10 µg/L. We investigated impacts of the FAR on drinking water violations temporally and geographically using the Safe Drinking Water Information System. Violations exceeding the MCL and the population served by violating systems were analyzed across the conterminous US from 2006 (onset of FAR enforcement) to 2017. The percentage of public water system violations declined from 1.3% in 2008 to 0.55% in 2017 (p < 0.001, slope = -0.070), and the population served decreased by over 1 million (p < 0.001, slope = -106 886). Geographical analysis demonstrated higher mean violations and populations served in certain counties rather than evenly distributed across states. The decline in violations is likely due to the adoption of documented and undocumented treatment methods and possibly from reduced environmental releases. Considering other studies that have shown decreased urinary arsenic levels in the population served by public water systems since the new standard, it may be inferred that the FAR is facilitating the reduction of arsenic exposure in the US.

Revising the index of watershed integrity national maps.

Watersheds provide a range of services valued by society, incorporating biotic and abiotic functions within their boundaries. Recently, an operational definition of watershed integrity was applied and indices of watershed integrity (IWI) and catchment integrity (ICI) were developed and mapped for the conterminous United States. However, these indices were originally derived using equally-weighted first-order approximations of relationships between anthropogenic stressors (obtained from the U.S. EPA's StreamCat dataset) and six watershed functions. In addition, the original calculations of the IWI and ICI did not standardize metrics across these differing scales, resulting in IWI and ICI values that are not directly comparable. We provide an example of how to iteratively update the stressor-watershed function relationships using random forest models and a nationwide response metric representative of one of the six watershed functions. Specifically, we focused on the chemical regulation function (CHEM) of IWI and ICI by relating a composite metric of chemical water quality from 1914 samples to land use metrics explicit to CHEM to refine the nature of these relationships (e.g., non-linear versus linear). The rate of nitrogen fertilizer, agricultural land use, and urban land use were found to be the three most important stressors predicting the national water quality response metric. Revision of CHEM values improved the prediction of several regional- to national-scale water quality indicators. In all cases, exponential decay curves replaced the original negative linear relationship for CHEM. Therefore, the original IWI and ICI values are probably over-estimates of the actual integrity of the Nation's watersheds and catchments. With these revisions, we provide updated national maps of IWI and ICI. The methods outlined here can be implemented iteratively as more and better data become available for all six of the watershed functions to elevate the accuracy and applicability of these indices to various land management issues.